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We study the linear conductance of single electron devices showing Coulomb blockade phenom- 
ena. Our approach is based on a formally exact path integral representation describing electron 
tunneling nonperturbatively. The electromagnetic environment of the device is treated in terms of 
the Caldeira-Leggett model. We obtain the linear conductance from the Kubo formula leading to 
a formally exact expression which is evaluated in the semiclassical limit. Specifically we consider 
three models. First, the influence of an electromagnetic environment of arbitrary impedance on a 
single tunnel junction is studied focusing on the limits of large tunneling conductance and high to 
moderately low temperatures. The predictions are compared with recent experimental data. Sec- 
ond, the conductance of an array of A'^ tunnel junctions is determined in dependence on the length 
A'^ of the array and the environmental impedance. Finally, we consider a single electron transistor 
and compare our results for large tunneling conductance with experimental findings. 



I. INTRODUCTION 

Tunneling of electrons in nanostructures is strongly 
affected by Coulomb repulsion. In systems containing 
metallic tunnel junctions the interaction can be described 
by the charging energy [nl Ec = e^/2C expressed in 
terms of a geometrical capacitance C. For weak tun- 
neling and temperatures well below Ec/ks, tunneling is 
suppressed by the Coulomb blockade effect. This regime 
is well explored experimentally P-pl , and the phenomena 
observed can be explained theoretically P-M by means 
of perturbation theory in the tunneling strength which is 
characterized by the classical high temperature tunneling 
conductance Gt- This approach breaks down for con- 
ductances comparable to or even larger than the conduc- 
tance quantum Gk = e^//i. When using these devices, 
e.g. as highly sensitive electrometers [O, in detectors 
p2| , or for thermometry fl^ , a large current signal is de- 
sirable meaning large tunneling conductance. However, 
higher order processes such as cotunneling 1 14 IJ] lead to 
a smearing of Coulomb blockade phenomena and a com- 
promise must be found in practice. While the strong tun- 
neling regime has been explored extensively by recent ex- 
periments |16|-|l9||, theoretical predictions remain limited. 
The theoretical work roughly splits into two groups. On 
the one hand, higher order perturbative results |20-2J] 
were successful in explaining some of the recent exper- 
imental data, yet, the analysis typically is restricted to 
conductances at most of order Gk- Based on the di- 
agrammatic expansion, partial resummation techniques 
were used to obtain nonperturbative results pq-|2q| , how- 
ever, for a restricted set of charge states. The arbitrary 
cutoff necessary in these latter theories limits their use 
for direct comparison with experimental findings. Fur- 
ther progress can be made by using perturbative renor- 
malization group techniques p9|,p0[. Apart from these 



approaches based on diagrams generated by treating tun- 
neling as a perturbation, a formally exact path integral 
expression [ pl[ including all orders in the tunneling con- 
ductance may serve as a starting point for analytical 
predictions |32-Bq] and numerical calculations [p7|-p9[. 
While perturbation theory in the tunneling term usually 
starts from states with definite electric charge, this latter 
approach employs the canonically conjugate phase vari- 
able and thus is well adapted to situations where the 
charge is smeared by thermal or quantum fluctuations. 
In this work we use the path integral approach to deter- 
mine the linear conductance of single electron devices in 
the semiclassical limit. Some limiting cases of the results 
presented here were published in short form previously 
|40|-[42[. Here we give a fuller account of the approach 
and apply it to a larger variety of systems. 

The paper is organized as follows: In Sec. II we intro- 
duce the Hamiltonians of a tunnel junction and of the 
electromagnetic environment, respectively. We then ex- 
plain the general method of calculating the linear conduc- 
tance from the Kubo formula with the help of a generat- 
ing functional. In Sec. Ill the case of a single tunnel junc- 
tion embedded in an electromagnetic environment of ar- 
bitrary impedance is considered. We use this example to 
derive the effective action characterizing the generating 
functional which is employed also in subsequent sections 
with adequate generalizations. We evaluate the conduc- 
tance in the semiclassical limit appropriate for high tem- 
peratures and/or large tunneling conductance and com- 
pare the results with experimental findings by Joyez et al. 
]l7t and by Farhangfar et al. ]18] . As a first extension of 
the method, we consider in Sec. IV a linear array of tun- 
nel junctions embedded in an electromagnetic environ- 
ment. The conductance of the array is determined in the 
high temperature limit. Specifically, we study the effect 
of the environmental impedance on the conductance and 



show that with increasing length of the array the influ- 
ence of the environment is strongly suppressed. In Sec. V 
we turn to a single electron transistor (SET). Here, we go 
beyond leading order in the semiclassical expansion and 
determine the conductance in dependence on the gate 
voltage. The findings are compared with experimental 
data by Joyez et al. nm for SETs in the strong tunneling 
regime. We conclude and discuss possible extensions in 
Sec. VI. 



II. MODEL AND GENERAL METHOD 

In this section we introduce the Hamiltonian for a sin- 
gle tunnel junction and model the electromagnetic en- 
vironment in terms of a set of LC circuits. A metal - 
oxide layer - metal tunnel junction consists of two metal- 
lic leads separated by a thin oxide layer [^Q . Provided 
the screening length in the metal is small compared to 
typical electrode and oxide barrier dimensions, one may 
introduce a geometrical capacitance C. The energy shift 
for an electron tunneling from one lead to the other is 
determined by the charging energy Ec = e^/2C. The 
corresponding Coulomb Hamiltonian reads 



HciQ) = ^, 



(1) 



where Q is the charge operator on the capacitance. The 
leads are described in second quantization by 
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where the tpa are quasiparticle energies, and a^^,^ and apa 
are creation and annihilation operators for states on the 
two electrodes, respectively. The indices p = k, q are 
longitudinal wave numbers and a is the channel index 
including transversal and spin quantum numbers. Pro- 
vided the tunneling amplitudes are small, we may de- 
scribe barrier transmission by a tunneling Hamiltonian 

0a 



Ht{(P) = y2 i'^kqaal^aqaA + H.C 



kqa 



(3) 



preserving the channel index <t. Here t^qa is the tun- 
neling amplitude and A the charge shift operator obey- 
ing A^QA = Q + e. Defining a conjugate phase Lp by 
[Q, ip] = ie, we may write 



A = exp{—iip). 



(4) 



The total Hamiltonian of a tunnel junction then reads 



Hj{Q, ^) = HciQ) + ^qp + Ht{^), 



(5) 



where the dependence on the charge and conjugate phase 
operators is made explicit to emphasize the similarity 



between the charging energy and a kinetic energy and 
between the tunneling Hamiltonian and an effective po- 
tential energy. 

The electromagnetic environment can be described by 
a Caldeira-Leggett model |45 
LC circuits 



as a linear combination of 
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coupled to the phase operator ip of the device. The pa- 
rameters of the LC-circuits are related to the environ- 
mental admittance by 



N 
Y{uj) = ^^-j— [S{^ + ^n) + S{uJ - UJn) 



(7) 



where the w„ = l/\/LnCn are the eigenfrequencies of 
the oscillators. A single electron tunneling device con- 
sists of tunnel junctions, capacitances, and admittances. 
The Hamiltonian of this system can be constituted from 
the elements discussed above. Then, the bosonic degrees 
of freedom of the admittance and the fermionic degrees 
of freedom of the electrodes may be traced out. In the 
next section we exemplarily derive the Hamiltonian and 
the effective action for a tunnel junction embedded in an 
environment of arbitrary impedance. 

Specifically, in this article we determine the linear con- 
ductance G = dI/dV\v=o of single electron devices. 
Here / is the measured current and V the applied volt- 
age. Within linear response theory one may use the Kubo 
formula for the linear conductance 



G{uj) = —- lim 
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dT e"'-^l'^^\T)I^^\Q)) , (8) 



where l'^^-' is the measured current and I^^-* a current op- 
erator determined by the coupling to the applied voltage 
V , see below. The v^ = 2TTn/h(3 are Matsubara frequen- 
cies. Correlation functions can be written as variational 
derivatives 



(/W(r)/(^)(r')) = 
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depending on auxiliary fields ^i. Here, the Hamiltonian 
H describes the system at vanishing external voltage, 
V — 0, and Tr is the Matsubara time ordering oper- 
ator. In subsequent sections we apply this generating 
functional approach to derive an explicit expression for 
the linear conductance in the semiclassical regime. 



III. TUNNEL JUNCTION WITH 
ENVIRONMENT 

A. Generating Functional 

We consider a tunnel junction characterized by the ge- 
ometrical capacitance C and the tunneling conductance 
Gt embedded in an electromagnetic environment. Via 
network transformations it is always possible to trans- 
form the environmental degrees of freedom into an ad- 
mittance Y(uj) in series with the junction biased by a 
voltage source V, cf. Fig. || In this subsection we obtain 
the effective action characterizing the generating func- 
tional introduced above. Readers familiar with path in- 
tegral techniques for single electron devices may directly 
proceed to the next subsection. 
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FIG. 1. Circuit diagram of a tunnel junction in series with 
an admittance. 

The Hamiltonian is given by _ff = Hj{Qj,ipj) + 
Hcmi^cui)- Here the phases (pj and ipcm are related 
to the voltages Vj and T4m across the tunnel junction 
and the admittance, respectively, by (pj = y-Vj and 
'fern — fKjm- Further, one has to take care of constraints 
for the variables imposed by the circuit. Using Kirch- 
hoff's law for the voltages, we find that the sum of the 
phases in the circuit loop in Fig. yj has to be constant, 
i.e. if J -\- Lpcm + i^ = const. ^ where we have described the 
voltage source in terms of an additional phase |g|] 



V'(i) = ^ / dt'v{t') 



(11) 



Similar relations hold for each loop of more complicated 
circuits. For an adequate handling of these constraints 
we start from the Lagrangian description, C = T — U . 
In general, the kinetic energy T is given by the sum 
of Coulomb energy terms and the effective potentials 
are the tunneling and environmental Hamiltonians. The 
constraints are naturally implemented by expressing the 
variables through generalized coordinates. Defining gen- 
eralized momenta in the standard way, one can derive the 
Hamiltonian via a Legendre transformation. To define 
conjugate momenta non-ambiguously, we use the discrete 
Caldeira-Leggett model and perform the continuum limit 



only afterwards. Shunt capacitors need to be treated sep- 
arately and will be discussed in Sec. V. Since ^{t) is con- 
trolled externally, the phase (^om rnay be eliminated in 
favor of ipj = (/3 and we may write 



^Je(Q,<^) =Hj(0,<^)+i?om(<^ + ^), 



(12) 



where Q = j^dC/dtp is the momentum canonically con- 
jugate to Lp. In the second term, we have absorbed the 
minus sign in front oi i.p + tp into the arbitrary definition 
of the sign of the phase of the environment. The current 
may be defined as the time derivative of the charge 



Q = t[H}E:Q]=It + I. 
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where 



It = ■^[Ht{(p), 
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is the current through the tunnel junction and 
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(15) 



the current through the admittance at vanishing exter- 
nal voltage. To determine the linear conductance (H), we 
choose the measured current /^^^ to be the current Jem, 
and l'-'^' follows from the coupling to the phase -0 in linear 
approximation: iJjE(Q, (^s, (Si/') = i?jE(Q,</') + ^I^'^^Sip. 
Via a unitary transformation U = exp(— zKii/^Q/e), it is 
always possible to write the external voltage partly as a 
shift of the phase variable ip. Using U'(pU = ip — Kitp 
and the general relation H' = U^HU + ihU^ -^U, we get 

Hje{Qi V) = -ffj(<9 + KlVC, ip - KiV^) + iJom(<y5 + K2V'), 

(16) 

where ki is an arbitrary shift and K2 = 1 — ki. Here 
we choose ki = so that the voltage couples solely 
to the environmental degrees of freedom and then get 
HcYa{p+5i}) — Hem{ip) + ^IemSip- Hcncc, in this case /^^^ 
coincides with the measured current /(^' — I^^^^ = I^m- 
To derive the path integral representation of the gener- 
ating functional (|lO|), we define 

~ e 

Hc-m{^)=Hcn,{p) - i{T)Ic-m{ip)^Hc-mVp - T^{'^)] + ind. 

n 

(17) 

where ind. denotes a (^-independent term that may be 
omitted. Further, we separate the exponential in Eq. (Hfl) 
into a free part Aq ihfi) for vanishing tunneling and a tun- 
neling part Axihp) according to 
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(18) 



where 
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= T.^-iIo''^'"°'-^'^ 



(19) 



describes the system in presence of the unperturbed 
Hamiltonian Hq ~ iJcm + He + Hqp. Using the series 
expansion of At{t) in powers of Ht and separating the 
trace in Eq. (pl]| ) into partial traces over the charge de- 
grees of freedom of the device, the quasiparticle compo- 
nents, and the environmental degrees of freedom, we ob- 
tain an expression of the generating functional as a sum of 
averages of the unperturbed system, cf. Q . Due to the 
Coulomb interaction Hamiltonian He in the unperturbed 
Hamiltonian Hq, contributions of a given order in Ht 
cannot simply be evaluated with the help of Wick's the- 
orem, however, the partial traces over the quasiparticle 
components are averages weighted with the free fermionic 
density matrix '^ eKp(—f3Hqp) and, accordingly, products 
of quasi-particle creation and annihilation operators in 
the interaction picture decompose into products of two- 
pair correlators. In the limit of large channel number 
Af = J2<7 1^1, only specific combinations of contrac- 
tions contribute that may be written in terms of two-time 
correlators of the tunneling Hamiltonian 
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with a real averaged tunneling matrix element t = t^qa- 



Here 



denotes the thermal average over the quasi- 



particles with Hamiltonian -ffqp. The time dependence 
in the interaction picture reads 

Ht{t) = exp (JiJqp) i^^Texp (- Ji^qp) (21) 

and 



K^{t)=Aq{-t)K^Aq{t). (22) 



Further, we have introduced the notation a"*" = a^, 
a~ = a, and A = exp(=pz(^). Performing the contin- 
uum limit for the longitudinal quantum numbers k and 
g, we find 

G{Ty) = \^a{r - t') [At(r)A(r') + At(r')A(r)] 



flGK 



(23) 



where Gt/Gk — '^ir'^i^M pp' is the classical dimension- 
less tunneling conductance with the densities of states 



p and p' at the Fermi level in the left and right elec- 
trode, respectively. In our approach the limit of strong 
tunneling is defined by A/" ^ 1, t^pp' ^ 1 such that 
Air^i^Mpp' ^ 1. Since for lithographically fabricated 
metallic tunnel junctions typically A/'-lO^, Gt/Gk can 
become very large, although each single channel is weakly 
transmitting only. 

The quasiparticle excitations generated by Ht are de- 
scribed by an electron-hole pair Green function ]20[| 
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(24) 



where the electron and hole propagate on different elec- 
trodes. D is the electronic bandwidth which may be 
set to infinity at the end of the calculation since D ^ 
EcksT. Due to analytic properties of thermal Green 
functions, wc may write 



a{T) 



1 
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with Fourier coefficients 
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Here and in the remainder of the article the absolute 
value is defined by 



z Re(z) > 

-z Re(z) < 



(27) 



which leads to a unique analytical continuation |}46| of 
the Fourier coefficients (P6[). Along these lines the partial 
traces over the quasiparticle components may be evalu- 
ated in terms of the tunneling kernel a{T). 

To proceed we need to consider next the partial trace 
over the charge degrees of freedom. It is convenient, 
to change to the phase representation and insert iden- 
tity operators J diy9r|<Pr)(<Pr| at each imaginary time 
shce r„ = ^n, n = . . . , iV with iV ^ oo. The 
charge shift operators in the interaction picture then 
become A (r) = exp(=p«(/jT-). Dividing the generating 
functional (|lC|) by the quasiparticle partition function 
trqp exp(— /3iJqp) which has no effect on the correlator 
d), we get 

^JeK]= I'D[^]\{ /'p[0„]cxp|-i5o[^,0„,C] 
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where S'o[<p, (/in,^] = Self] + Sem['^,4>n,i] contains the 
environmental and the Coulomb actions specified below. 
Since the integrand is invariant under exchange of an ar- 
bitrary pair of variables we may extend the integrations 



to L dTi {i — 1, . . . , 2m) and compensate the larger in- 
tegration region by a factor l/{2m)\. Further the sum 
over pairs leads to a factor (2m — 1)!!. Interchanging 
integrals and product we get 






(29) 



where the effective Euclidean action splits into three 
parts 

SjEif, </»», C] = Self] + StM + Scmi'P, <f>n, £,] ■ (30) 

Here 
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describes Coulomb charging and 
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quasi-particle tunneling across the junction. The envi- 
ronmental action is given by 
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The remaining trace over environmental degrees of free- 
dom in Eq. (p9|) can be evaluated exactly E^ leading to 
a quadratic nonlocal action 



Sy[^,^]^^J drj dT'k{T-T') 



V{r) ~ -£,{t) 



Mr') + -^ar') 



(34) 



where the kernel ^(t) can be written as a Fourier series 
(Eq) with coefficients 
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Here Y{s) is the Laplace transform of the environmental 
response function F(i), c/. Rcf. |47|. Due to causality, 
for Re(s) > 0, one may write Y{s) — Y{is) where 1^((jj) 
is the frequency dependent admittance (H) of the envi- 
ronment. 



This way the generating functional reads 

^jEK]=y2?[^]exp|-i5jE[^,C]| 
with the effective action 

'S'jE [<^, C] = Sc [f] + St W\ + Sy[^,^]- 
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The explicit form of the generating functional serves as 
a starting point to calculate the correlator in the next 
subsection. 



B. Conductance 

We now perform the functional derivatives in Eq. (0) 
explicitly and get for the correlator p3 

(/e„.(T)/em(0)) = ^ I I?M exp | - ^^JE [^, 0] | 



X ( 2 — fc(T) +/„„[¥>, T]/e,n[(^,0]j , (38) 

where Zje — ^je[0] denotes the partition function. The 
current functional Iem[fiT] arising as variational deriva- 
tive of the effective action (|37|) reads 



2p r"-^ 
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The conductance (ra) now splits into two pieces. 



GjE(u;) = GW(c.)-fGfE)(c.) 
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to the second term in Eq. (pq). The cxphcit form of 
the auxiliary functional i^[(/3, j/„] in terms of the Fourier 
components ip{vrn) reads 
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m— — oo 



So far no approximations have been made and Eqs. (|40|) — 
(H) give a formally exact representation of the linear 
conductance. To proceed we evaluate the path integral 
(Es) in the semiclassical limit. 



C. Semiclassical Limit 

The classical trajectory of the phase (p is defined by 
SSjE['f,0]/S'P — 0, and we obtain from Eq. (pTj ) (p — 
ip^ = const. Since the action is invariant under a global 
phase shift, we may put (p" = 0. Writing the action in 
terms of Fourier coefficients of the phase 
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and expanding in powers of (/?(^'^), we get 
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is the second order variational action with the eigenvalues 
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Here 



Goiiyn) = Wn\C + Gt 



(49) 



describes the tunnel junction as a capacitance in paral- 
lel with an Ohmic resistor characterized by the classical 
tunneling conductance. Further 
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is the variational action of order 2k. The summation over 
the Hi is over all integers with rii = omitted. Neglecting 
sixth and higher order terms, we get from Eq. (p3|) 
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The convergence of this expansion depends crucially on 
the eigenvalues (E3) . To estimate the range of validity of 
the truncated series, we write the smallest eigenvalue in 
more appropriate units as 



Aje(;^i) 



27r^ Gt + Y{vi) 
PEc Gk 



(52) 



This eigenvalue has to be large compared to 1, and 
we see that the expansion is useful for large conduc- 
tances Gt + ^^(1^1) ^ Gk and/or high temperatures 
l3Ec ^ Stt'^. Hence, we effectively expand in powers 
of 



e = Min 
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Performing the limit w„ — > uj 
tinned eigenvalue (Es) reads 



AjE(-«tj) = -iuj^- [Go{uj) + Y{uj)] , 



where 



Go{lu) — Gt — i^G 



(54) 



(55) 



is the analytic continuation of the Laplace transform of 
Eq. (p9|). The relative minus sign of the capacitive term 
is due to the usual definition of the Fourier transform in 
quantum mechanics, the electro-technical convention is 
obtained by replacing oj — > — w. For small frequencies the 
analytically continued eigenvalue (B3) is no longer large 
compared to 1 and we are faced with a problem of order 
reduction. In the limit ivn -^10 + 18 each l/AjE(^'ri) in 
Eq. ( pl|) becomes of order 1 while the l/AjE(t'm) factors 
for m ^ n are not analytically continued and remain of 
order e. The correction term of order e^ in Eq. (|5l|), that 
is the term proportional to GtJGk-, becomes of order e 
after analytic continuation. Hence we loose one factor 
of e. Generally, one finds that the higher order varia- 
tional actions ( pOf ) include at most one l/AjE(^'n) factor 
and consequently are reduced at most by one order in e. 
Thus, products of the form 



s'An'p]s'Anv']---s': 
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(56) 



of quartic or higher-order variational actions, as they 
arise from an expansion in powers of ip[vk), give quan- 
tum corrections of order ^'■^1+^2+ ■■■+ki g^j^^ after analytical 
continuation of order ^^^1+^2+ ■■■+ki-i ^^-^^ of higher orders. 
This proves that the terms of the expansion of i^jE given 
explicitely in Eq. ( pl| ) suffice to calculate the leading or- 
der quantum corrections. After performing the analytical 
continuation we get 
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with the quantum correction factor 
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Hence, for the total conductance we may write 
GjE(a;) 
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(59) 



G,a{oj) + Y{Lo) 

with an effective linear conductance of the junction 

Geff (w) = Gt [1 - U{uj)] - iLuC. (60) 

This describes a linear element G*{uj) — Gt[1 —U{uj)], 
depending on the whole circuit, in parallel with the geo- 
metrical junction capacitance C, cf. Fig. ga. The general 
form ( p9| ) is valid only to first order in e. A systematic 
treatment of higher order contributions does not allow 
for a description of the tunnel junction in terms of an 
effective linear element. However, a partial resumma- 
tion of higher order terms according to a self-consistent 
harmonic approximation p? 48 1 leads again to the form 

(ID- 
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FIG. 2. Effective circuit diagrams for a tunnel junction in 
the semiclassical limit a) for arbitrary frequency and b) in the 
low frequency limit. 



D. Results and Comparison with Experimental Data 

For further discussion and comparison with experi- 
mental data we restrict ourselves to ohmic dissipation 
Y{ll!) ~ Y. The effective linear element ( pOJ ) then reads 



G^M 
Gt 



1 



■ijj{l + u + Lu) -il;{l+uj) 



■0(1 +U + LU) --0(1 +U) 



fiEc 



(61) 



where '0 is the logarithmic derivative of the gamma func- 
tion and 



u = g- 



PEc 
27r2 ' 



h(5 
2ttI'' 



(62) 



are auxiliary quantities. We also have introduced the 
dimensionless parallel conductance 



g^{GT + Y)/G 



K- 



(63) 



The quantum corrections depend only on this combina- 
tion of conductances. The real and imaginary parts of 
G*{uj)/Gt are depicted in Fig. || for j3Ec — 1 and vari- 
ous values of g. 

The quantum corrections are most pronounced at zero 
frequency and disappear nonalgebraically for large uj 
and/or u, due to the logarithmic behavior of the psi- 
function for large arguments. 
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FIG. 3. Real and imaginary parts of G*{lo)/Gt in the 
ohmic damping case for PEc = 1 and various values of the 
dimensionless conductance g in dependence on the dimension- 
less frequency Q. = fiuj/2iiEc- 



For the dc conductance we get from ( pi] ) 



G*{uj = 0) 
Gt 



7 + -0(1 + u) 



ij'{l + u) 



l3Ec 



(64) 



which coincides with our previous result Effl. In partic- 
ular, in the limit of a very low resistance environment, 
the total conductance (|5£ ) approaches the classical limit 
nonanalytically, cf. FigTE, leading to the asymptotic ex- 
pansion Eg] 

G.M^ = 0) 



Gi 



1 



2^ In 
Y 



Gk 
Y 



O 



HPEc)\G 



K 



Y 



(65) 



On the other hand for moderate to large environmental 
resistance, we may expand Eq. (p4) with respect to u 
leading to a total conductance 



Gje(c^ = 0) 
GtY 



Gi 



Y 



1- 



y PEc 

G + Y 3 



^0[i(3Ecf,uPEc]\. (66) 



This approximation correspond to the dotted line in 
Fig. y and remains analytic in the limit of large envi- 
ronmental conductance where it obviously fails. 




Gk/Y 

FIG. 4. The ratio of the total dc conductance Gje(O) and 
the classical conductance Gci ~ GtY/{Gt + Y) shown vs. the 
dimensionless environmental resistance Gk/Y for Gt ~ Gk 
and (3Ec = 1/4. The solid line is the result (^) for oj = 
and ohmic damping, and the dotted line corresponds to the 
approximation (pq) for moderate to small Y. 

In Fig. ^ we compare our prediction (|6j) with recent 
experimental data by Joyez et al. [|l7| for dimension- 
less conductance g = 4.2 and 23.8 (upper plot) and by 
Farhangfar et al. |Q for g = 4.52 and g = 34.2 (lower 
plot). Fig. ^ shows that in the limit of large conductance 
we are able to explain the whole range of temperatures 
explored experimentally, whereas for moderate conduc- 
tance only the high temperature part is covered by the 
semiclassical theory. Here, the parameters g and Ec have 
not been adjusted to improve the fit but coincide with 
those given in the experimental papers. 
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FIG. 5. The renormalized conductance (p4) versus the di- 
mensionless temperature compared with experimental data 
by Joyez et al. wA for dimensionless parallel conductance 
g — 4.2 and 23.8 (upper plot) and with experimental data 
by Farhangfar et al. mM for g — 4.52 and 34.2 (lower plot). 

We conclude this section with some remarks on the fre- 
quency dependence of the conductance that has not been 
studied experimentally, so far. For snrall frequencies we 
may expand the result (|6l|) and write 



G*{lu) = G*{u; = 0)- iujC* + 0(cj2), 



(67) 



where G* leads to a renormalization of the junction 
capacitance C. The renormalized capacitance GcS — 
G + C* reads 



Gcs 
G 



Gj' 
Gk 



2V''(l + w) 



i:"{l + u) 



{PEc 



47r4 



(68) 



The correction shows a quadratic dependence on PEc 
and therefore is suppressed at high temperatures. It also 
vanishes linearly for large conductance g due to the an- 
alytical properties of the psi function. The semiclassi- 
cal treatment covers only the region of weak Coulomb 
blockade. Whereas for small tunneling conductance low 
temperatures imply strong Coulomb blockade, these ef- 
fects are suppressed for highly conducting tunnel junc- 
tions and the semiclassical theory is restored. A closer 
examination of Eq. ( |6l|) shows that for g '3> 2\n{f3Ec) 
the quantum corrections are always small. For fixed 



g ^ 1, our predictions are therefore valid for a very 
large range of temperatures covering in fact the en- 
tire range of parameters presently attainable experimen- 
tally for metallic junctions with strong tunneling p7|-fi9[. 
Fig. @ depicts the real and imaginary parts of G*{lu)/Gt 
for g = 60 and various temperatures. With decreas- 
ing temperature the real part shows for a; = a loga- 
rithmic decrease, G*/Gt = 1 — 2 ln(/3 i?c)/5j as long as 
ksT ^ Eccxp{—g/2). Thus, for large conductance the 
semiclassical treatment is an effective high temperature 
expansion valid for ksT large compared with the renor- 
mahzed charging energy E^ w i^c exp(— g/2) p3| , p6[ . 
Moreover, the analytical form of the quantum corrections 
indicates that Coulomb blockade survives for arbitrary 
large conductance but becomes strong only for tempera- 
tures below E^/ks- 

In the limit T — j 0, g — > cxd such that /csT 3> 
Ec exp{—g/2), the imaginary part of G*{u!) becomes a 
step function of width 2TT/g, cf. Fig. ^, leading to a di- 
vergent renormalized capacitance of the form Ccs/C ~ 
PEcGt/Q{Gt + Y). The linear dependence of Ccg on 
/3 starts already at very high temperatures, cf. Fig. |j, 
and only saturates for (3Ec of order exp{g/2). The large 
renormalized capacitance is a strong tunneling effect due 
to multiple electron tunneling and is found likewise for 
non-Ohmic environmental impedances. In Fig. ^^we show 
the renormalized capacitance Ces/C for Gt/Gk — 20 
and various values of Y/Gk in dependence on the dimen- 
sionless inverse temperature l3Ec- Note that the linear 
behavior of the capacitance starts already near (3Ec — 1 • 
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FIG. 6. Real and imaginary parts of G*{lo)/Gt in the 
Ohmic damping case for dimensionless conductance g — 60 
and inverse temperatures fiEc = 20, 40, 80, and 160 in de- 
pendence on the dimensionless frequency Q — uTi/2-nEc- 

The renormalized capacitance describes the frequency 
dependence of the conductance for small frequencies 
ujCeS <C ■nGTJg., cf. Fig. |[ Rewriting this inequality 



we get u! <C QkBT/U « IQ-'^-'^THz, where T is the temper- 
ature measured in Kelvin. Thus for all accessible tem- 
peratures the frequency range of strong 1// noise can be 
avoided, and the effect predicted should be clearly ob- 
servable experimentally. 
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FIG. 7. Renormalized capacitance Ces/C in the Ohmic 
damping case for tunneling conductance Gt/Gk = 20 and 
various environmental conductances Y/Gk = 1, 5, 10, and 20 
in dependence on the dimensionless inverse temperature jSEc . 



IV. ARRAY OF JUNCTIONS WITH 
ENVIRONMENT 

A. Generating Functional and Conductance 

As a first extension of the method, we now consider lin- 
ear arrays of A'^ tunnel junctions embedded in an electro- 
magnetic environment. The junctions are characterized 
by classical tunneling conductances Gj and geometrical 
capacitances Cj in parallel. Like in the previous section, 
the environment can be transformed into an admittance 
y(w) in series with an array of junctions biased by a 
voltage source V, cf. Fig g. We start with a Lagrangian 
description depending on phase variables (pj of each junc- 
tion j = 1 . . . A^ and an environmental phase </Jcm with 
the constraint X^i^i "^i + V'ein + V' = const., where %/: 



describes the applied voltage and is given by Eq. (|1] 
Using the ipj, j = 1 . . . N as generalized variables we find 
for the total Hamiltonian 



N 



N 



^Ae({Qj}, Wj})=Y,H3{Qj,Vj) + i/em ^ P, 



i=i 



3=1 



(69) 



with the junction and environmental Hamiltonians de- 
fined by Eqs. (j!]) — (^)- We follow the analysis in the 
previous section and first derive a formally exact expres- 
sion for the linear conductance. As measured current P^' 
we choose again the current flowing through the environ- 
mental impedance given by Eq. (|l5|) with ip replaced by 
J2j=i'Pj and Q by J2j=iQj^ respectively. The second 
current operator /'^' is determined by the linear cou- 
pling to V and we get /^^^ = I^^^ = /em- Following the 



lines of reasoning in the previous sections, the generating 
functional is found to read 



^aeK] = j D[{^,}] exp |-i 5AE[{</5,},e]} , 



(70) 



with the effective Euclidean action 



Sa'e.[W]},S,] = Sy 






N 



E=-'^^'^ +E^^[^^]- (71) 



Here Sy was introduced in Eq. (|34j) and Sj[lpj\ — 
^fVPj] + ^JiVj] describes the j'th junction where the 
Coulomb action S'^ and the tunneling action Sj are 
given by Eqs. ( |3l| ) and (|3^, with the replacements Gt -^ 
Gj and C ^ Gj. Performing the functional derivatives 
explicitly, the current-current correlator is found to be of 
the form ( p8[ ) with the replacement Zje —^ Zae = Zae [0] 
and the appropriate action S'AEll'ySi}, 0]. Further, the 
current functionals Iem[^,T] now depend on the sum of 
phases, ip — > X)j=i¥'i' ^^'^ the functional integral is 
defined over all configurations of the phases ipj. As in 
Eq. (Bl|) the first term can be handled exactly, and we 

— (2) (2) 

find Gae(w) = Y{uj) + G)^g(w), where G)^g(a;) given by 
Eq. (H) and Eq. (||) with F[(^,z/„] -> i^Ef=i '/'j, '^n]- 
So far no approximations have been made and the result 
follows from a straightforward extension of our findings 
for a single junction. The qualitative difference lies in 
the topological structure of the phase configuration space 
and becomes clear when one evaluates the path integral. 
Again, we restrict ourselves to the semiclassical limit. 
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FIG. 8. Circuit diagram of an array of N tunnel junctions 
in series with an admittance Y{lj). 



B. Semiclassical Limit 



torus. Accordingly, the path integral is over all configura- 
tions of the phases ipj with ipj (hpi) = ipj (0) -|- 2TTkj where 
the winding numbers kj are integers. On the other hand, 
the environmental phase is extended and conjugate to a 
certain linear combination Q of the Qj . Since the envi- 
ronment transfers charges continuously, Q is not quan- 
tized and the path integral over the environmental phase 
is over all configurations with tpemiO) = V'ein(7i/3)- Rather 
than making the canonical transform to the charges qj 
{j — 1, . . . , N~l), Q and the conjugate phases explicitely, 
one finds that equivalently we may integrate over all con- 
figurations of the phases ipj with ipj {h(i) = ipj (0) -I- 2Trkj 
where the integer winding numbers obey the constraint 

At high temperatures the classical paths are straight 
line flips (p, (t) ~ (p^^ + v^^t running from if^, to 
(/j*^ 4- 2i:kj. The action is invariant under global shifts 
of the ipj and we may set ipR = for all j — 1, . . . ,N. 
All paths with winding number kj ^ are exponen- 
tially suppressed by the classical action contribution 
Sf « TT^k]/l3Ec^ + \kj\Gj/2GK- Thus, to obtain the 
leading order quantum corrections, we may restrict our- 
selves to winding numbers kj — for j — 1 . . . N. Finite 
winding numbers are considered in the next section where 
we focus on the single electron transistor and go beyond 
the leading order quantum correction. 

The action may be expanded in powers of the Fourier 
coefficients (pi{i'n) yielding a result of the form (Q) where 
the second order variational action reads 



sUWA] = s'y 



^N 



Lpj 



N 



E^°[^. 



(72) 



The environmental contribution is given by 

oo 

S'^[ip]^hJ2>^Y{'^n)m'^n)\^ (73) 

n=l 

with the eigenvalues 






(74) 



The second order variational tunneling action for junc- 
tion j reads 



Determining the classical path one has to take into 
account the topological structure of the configuration 
space. The phases ipj, j = I, . . . ,N are canonically con- 
jugate to the charges Qj on the junction capacitances. 
Now the array has A^ — 1 metallic islands in between 
the junctions carrying the charges Qj = Qj — Qj+i for 
j — 1, . . . ,N — 1. These island charges are quantized 
in units of the elementary charge e, and the phases ipj 
canonically conjugate to the qj are compact, i.e., the con- 
figuration space of the phases is a {N — l)-dimensional 



S°[ipj] = hJ2^ji^n) |<Pi(l^n)r 



with the eigenvalues 



n=l 



h^(3 



^ji^n) = ^5-|l'„|G°(l^„), 



(75) 



(76) 



where G^Avn) is given by Eq. (E3) adapted to a junc- 
tion with capacitance Gj in parallel with an Ohmic resis- 
tor 1/Gj. The higher order variational actions are given 



10 



by straightforward extensions of Eq. (|5^) . Expanding in 

powers of the higher order terms >5'^e[(,5], A: = 2, 3, . . ., we ^AEi^'nj 

are left with expectation values of products of the phase 

variables ^^ (;/„). It is now useful to define a Gaussian 

average 



Wo = 



1 






oo N 
ri=l j = l 



e^p (-^SUW,}]\ X 



(77) 



with the Gaussian partition function Z^g defined by the 
requirement (l)o — 1. The difference between Z^e ^^"^ 
the full partition function Zae is of order {(iEcY ^-nd 
may be neglected here. Due to the Gaussian form of 
the measure exp[— S'^g[{(/7j}]/?i], the averages of prod- 
ucts of Fourier coefficients ipj{i/n) decompose into sums 
over products of two point expectations. For two differ- 
ent phase variables I ^ I' we obtain 



{ipiiiyn)^l'{l^7n))o = -Sn,- 



N+1 

n '^^j^^") 



K^n), (78) 



with 



Af-l-lAf-l-l 



(79) 



Here and in the remainder we define XN+i{i^n) = Ay(i'n) 
and note that summation and multiplication indices run 
from 1 if not otherwise specified. For phase variables of 
the same junction we find 



{ipi(yn)Vl{^m))o 



1 



A«(t.„) 



where 



AW(z.„) = ^|i.„|G;(j^„) 



(80) 



(81) 



plays the role of an effective eigenvalue for phase fluctua- 
tions in junction I with all other phases ^j, j =^ I already 
traced out. Here, 



GK^„) = G?(^„) + 



N 



YiWnl) 






(82) 



may be considered as the Laplace transform of an effec- 
tive response function describing the circuit seen from 
junction I, i.e., a series of A'^ — 1 junctions and an envi- 
ronmental impedance in parallel to junction I where the 
junctions are described effectively by linear elements. In- 
cluding the fourth order variational derivative of the ac- 
tion, we get as a generalization of Eq. (O) 



^ "^ 1=1 ]^i 



Gi 2/3 
Gk A(t^„) 
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E 



a{vn) - a{vm) 



A(') (,.„,) 



(83) 



Now, the convergence of the expansion depends on the 
effective eigenvalues (p^). To estimate the range of va- 
lidity, we consider the smallest eigenvalues which at high 
temperatures are given by X^^^vi) « 2t:'^ / j3Eci. Again 
the analytic continuation gives rise to a reduction of the 
order of the quantum corrections in the expansion param- 
eter. For contributions with vanishing winding number 
the arguments given in the previous section apply like- 
wise to the present problem. On the other hand, for lower 
temperatures one has to take into account winding num- 
bers kj 7^ and finds that some of the eigenvalues tend 
to zero. The marginally stable fiuctuation modes lead to 
a breakdown of the simple semiclassical treatment. The 
appropriate extension of the semiclassical approximation 
was discussed elsewhere [Q. The topological structure 
of the phase space leading here to a breakdown of the 
simple semiclassical approximation at low temperatures 
even for large conductance is the main difference between 
the single tunnel junction with environment and circuits 
containing many junctions. As a result one finds that 
the truncated expression ( p3[ ) is valid up to first order 
in e = Max {PEc^ : j = 1, . . . ,N). After the analytical 
continuation i/„ — > —iuj + 6 we can write the total con- 
ductance in the compact form 

-1 -1 
1 x^ 1 



Gae(w) 



Y{lo) 
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E 



1 G^eff(^) 



(84) 



describing A^ -I- 1 linear elements in series: G^fj (a;) with 
j — 1, . . . ,N and the admittance Y{uj). Here, the G^g(w) 
are of the form (pQ) where the auxiliary functions Uj are 
given by Eq. (^^ with Aje replaced by A^^^ introduced 
in Eq. (pU). This is a straightforward extension of the 
result in the previous section valid to linear order in e for 
arbitrary conductances Gj and admittances Y(uj). 



C. Discussion of Results 

For a more explicit discussion of the results we con- 
sider now N identical junctions Gj — G and Gj — G. 
The eigenvalues (pH) then read 



Xil^n) ^ X^'H'^n) 



r/3. 



G%iy„ 



Y{W„\)G°Ai^„) 



(/v~i)y(|z.„|) + G"(i.„) 



(85) 
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and coincide for all junctions. For the total conductance 
of the array (|8J) we obtain 

G,siLu)Y{Lu) 



where 



GH = 



G,s{^) + NY{lu)' 



Gcff (w) = G[l - U{oj)] - iuC 



(86) 
(87) 



is the effective linear conductance of one junction. In this 
order each junction can be described by a linear element 
G*[ijj) = G[l -~U{u})\, depending on the circuit, in par- 
allel with the geometrical capacitance C as depicted in 
Fig. pK. To proceed we consider an Ohmic environment 
y(w) = Y and find for the effective linear element 



G*{uj) __^ ^{N -l)[i}j{l + UT+Co)-Tp{l + i:b)] ^ '^{l + UN + Co)-ilj{l+Co) 



G [^ ut ujv 

J 7r2iV ' 
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where 



UN 



G + NY (3Ec 



G 
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27r2 



UT 



G (3Ec 
Gk 27r2 ' 



h(3 

27ri' 



(89) 

are auxiliary quantities and Ec = e^/2C is the charging 
energy for one junction. For A^ = 1 we recover the results 
of Sec. Ill, of course. For a large array with A^ ^ 1, terms 
in Eq. ( p8| ) containing un drop out, and the quantum 
suppression becomes independent of y. Furthermore the 
high-temperature anomaly, c/. Fig. H, is now a 1/A^ effect 
and the limiting result for A^ — > oo is analytic. 
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FIG. 9. Renormalized conductance G* /G of an array 
of A'^ = 20 tunnel junctions in dependence of PEc for 
Y/Gk = 20 and various tunneling conductances G/Gk- 

For small frequencies the effective element behaves like 
an Ohmic resistor \/G*{uj = 0) with a renormalized ca- 
pacitance in parallel. The dc conductance is given by 

G*{lo = Q) 



G 

1-[{N 



7 -f -0(1 -f Ut) 



Ut 



i:'{l + UT) 
PEc 



7-t-V(l + ww) , ,//-, , X 

+ +1/^(1+ mat) 

UN I Att^ 

For N — 2 and Y/Gk -^ oo this reduces to 
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+ 0'(1 + mt) 



PEc/2 



(90) 



(91) 



coinciding with earlier findings for the symmetrical SET 
p5| , [ll| . On the other hand, in the limit of large A^ and 
moderate G~Gk, Eq- ( |90| ) reduces to 



G* 

IT 



N - 1 (3Ec 



(92) 



again in accordance with earlier findings MS] derived from 
rate theory for small tunneling conductances. To discuss 
the strong tunneling corrections, we show in Fig. M the 
renormalized conductance of an array of A'^ = 20 tun- 
nel junctions in dependence on (3Ec for Y/Gk — 20 and 
various tunneling conductances G/Gk- We find that the 
weakly conducting case G/Gk = 0-1 perfectly coincides 
with the limiting formula (92) (both depicted by the solid 
line) whereas for larger tunneling conductances strong 
deviations from this behavior appear. 

Experimentally one is interested in the dependence on 
the array length at fixed classical series conductance. In 
Fig. ^we show G* /G for various A^ whereby G increases 
with the array length to keep the total classical series 
conductance constant. Whereas for A^ < 5 the renor- 
malized conductance depends strongly on Y, it becomes 
independent for large A^. 
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FIG. 10. Renormalized conductance G* /G of an array of 
N = 1,2, 5, and 10 equivalent tunnel junctions leading to the 
same classical series conductance for jSEc = 1 as a function 
of the inverse environmental conductance Gk/Y. 
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The comparison with available experimental data | |5C| ] 
is complicated by the large number of parameters, in par- 
ticular, the charging energy differs from sample to sam- 
ple with different array length. To test our predictions 
at least in the perturbative regime, we compare with the 
results of a master equation approach [ pO[ based on the 
P{E) theory P|. In Fig. O we show the zero bias dip 
1 — G* /Gt in per cent for an array of length A^ = 20 and 
l3Ec = 0.0442 in the limit G ^ 0. We find good agree- 
ment between the numerical calculations by Farhangfar 
et al. Q and the analytical semiclassical result (^. 

The renormalized capacitance CeS includes the linear 
part in cu of G*{lu) and the geometrical capacitance C 
and reads 







x<iN-l) 



^-2V''(1 + ut) 



Ut 



>"(1 + mt) 



^'^^^T^-n^^^m 



(93) 



control voltage Ug shifting the electrostatic energy of the 
system continuously. The important energy scale is the 
charging energy Ec = e^/2C with the island capacitance 
C = Ci + C2 + Cg. For weak electron tunneling, Ec is 
the energy needed to charge the island with one excess 
electron at vanishing gate voltage Ug = 0. Due to the pe- 
riodicity of the Hamiltonian in Ug, the conductance is a 
periodic function with period 1 of the dimensionless gate 
voltage Ug = UgGg/e 0. 
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FIG. 12. Circuit diagram of the single electron transistor. 



showing a quadratic dependence on PEc ■ The renormal- 
ization is suppressed at high temperatures and also van- 
ishes linearly for large conductance in accordance with 
the behavior of a single tunnel junction with environ- 
ment. 




FIG. 11. Zero bias dip 1 — G* /Gt in per cent for an array 
of length A'^ = 20 and PEc = 0.0442 as a function of the in- 
verse environmental conductance Gk/Y in the perturbative 
limit compared with a numerical master equation approach 
by Farhangfar et al. [BOl. 



Following the lines of reasoning in the previous sections, 
we start with the Lagrangian description with phases 
(^1, (f2 across the tunnel junctions and Lpg across the gate 
capacitor. Here, we treat the shunt capacitor Gg ex- 
plicitely and do not introduce an effective environmental 
impedance. For the circuit depicted in Fig. O there are 
two independent circuit loops leading to the constraints 
<y2i ~ <<52 — "0 = const, and ipi — (pg — ipg — piip = const., 
where V is the phase (jl^) of the transport voltage V 
and Vg the corresponding phase of the gate voltage Ug. 
Eliminating ip2 and (pg in favor oi ip = ipi, we find 



-^SET 



HciQ + c 



-m 



, + V{G2 



0) + 



qp 



qp 



(94) 



The Coulomb Hamiltonian Hq is given by Eq. (|i|) and 
the other terms are defined by Eqs. (||) and (|^) with cor- 
responding labels. After a unitary transformation one 
finds equivalently 



i^SET = Hc{Q + eng) + H^r^\ip + mip) 

+h!^\p-k2^p) 



qp 



^(2) 

qp 



(95) 



V. SINGLE ELECTRON TRANSISTOR 



A. Generating Functional and Conductance 



with an arbitrary shift parameter ki and K2 = 1 — ki- 
Here we introduced a shifted dimensionless gate voltage 



ena 



UgGg + V{G2+PlGg-KlG). 



(96) 



The SET consists of two tunnel junctions with tun- 
nehng conductances Gi, G2 and capacitances Gi, G2, 
respectively, biased by a voltage source V, cf. Fig. |l^. 
The voltage may be split among the branches in piV 
and P2V with (pi + p2 = 1). The island in between the 
junctions is connected via a gate capacitance Gg to a 



In the sequel we restrict ourselves to zero frequency where 
the expectation values of the current operators through 
both junctions coincide (/i)(w = 0) = (-^2)(^ — 0). 
The first current operator /^^^ may then be chosen as 
an arbitrary linear combination J*^^^ = ei/i — £2/2 (with 
£1 + 62 = 1) of the tunneling current operators /i and 
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/2 through junctions 1 and 2, respectively. The relative 
minus sign conies from the opposite directions of Ii and 
I2, which are both positive for flux onto the island. The 
second current operator /'^' = kiIi — K2I2 is determined 
as above by the linear coupling term to the transport 
voltage V. The dependence of the Coulomb Hamiltonian 
on the transport voltage may be removed by a gate volt- 
age shift and thus need not be considered. Moreover, 
this coupling would lead to a displacement current con- 
tribution vanishing ai cu == 0. To evaluate the current- 
current correlator we employ the generating functional 
(p^. The current operators through the junctions are 
given by Eq. (|l4|), for j — 1,2. To derive the generating 
functional we define 



H, 



U) 



H^4\v)~iM^)- 



(97) 



The new tunneling Hamiltonians are of the form (0) 
where A is replaced by [1 + ie£,j{T)/h]Aj. With these 
replacements we get for the generating functional 



^SEtKi,6] = / D[ip]exp l--SsET[(P,^l,i2] 



(98) 



where the effective action reads 

SsetI^P, C] = ^s'etM + Sf[ip, a] + S^[ip, 6]- (99) 
The first term on the rhs 



-'SETL 



dT 



AEc 



+ ihngip{T) 



(100) 



describes Coulomb charging of the island in presence of 
an applied gate voltage. The effective tunneling action 



SJ[V,^A = ~^ 



hp 



dT 



np 



Gk Jo 



l + ^-^,{r') 



dT'a{T~T') l-i-a(T) 



^i[ip(r)-ip(T')] 



(101) 



describes quasi-particle tunneling through junction j 
with the kernel a{T) given by Eq. (|2q). For vanishing 
auxiliary field (,j = 0, the action reduces to the single 
electron box action |3^ S'set[</5] = <S'set[</', 0,0] = S'box 
depending only on the parallel conductance Gii — Gi + 
G2- One has 



^h0 



f^ pap f^np 

Sf[^,0] + S^[^,0]^2-^ dr dr' 

Gk Jo Jo 



a(r — r') sin 



(102) 



Thus, the generating functional at vanishing auxil- 
iary fields gives the box partition function ^set = 
■^SEt[0, 0] = .^box- Performing the variational derivatives 
explicitly, we get for the correlator 

{lAr)lAr')) = (/,(T)/,(r'))^<5„-- + {I,{r)I,,{r')f . 

(103) 



Since the auxiliary fields are in the argument of an ex- 
ponential, there are two contributions. The first term 
comes from the second order variational derivative of the 
action and reads 

a.(T)/,(r'))^-47rG,a(r-r')-^ /i?[^] 

^SET J 

expL^SsET[^]\ cos[^(r) - ^(r')]. (104) 



The second term in Eq. ( |103| ) involves a multiplication of 
two current functionals arising as first order variational 
derivatives of the action 

'^K -^SET J 

exp|-i5sETM|/[^,T]/b,r'], (105) 
with the current functional 



2e 



np 



I[ip,T] = — dT'a{T-T')sm[(p{T)-ip{T')]. (106) 



Taking into account that {Ij{t)Ij{t'))^ /Gj and 
{Ij{t)Ij'{t'))^ /GjGj' depend only on the parallel con- 
ductance G|| — G1+G2 and thus are independent of the 
indices j and j', the conductance may be written as 

G = eiKi{GiE + GlF) - {eiK2 + e2Ki)GiG2F 

+e2K2{G2E + GlF), (107) 

where 

1 f'"^ (/i(r)/i(0))-^ 

E= lim — Im lim ' -^--''^'^^ ^^ > ^^^^^ 



Gi 



(108) 



and 



F = lim - — Im lim 



+0 huj 



iui — fu+id J Q 



np 



dr ( 



,r{li{r)h{0)y 



Gi 



(109) 



Since the conductance does not depend on the specific 
choice of the parameters ej and kj , we then find that 



GsET — GciE, 

with the classical series conductance 

G1G2 



Gr] — 



G1 + G2 



(110) 



(111) 



This is a formally exact expression for the linear dc con- 
ductance. To proceed, we make explicit the sum over 
winding numbers k of the phase and write the correlator 
(104) in the form 
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ip{hl3)=2Trk 



(/i(t)/i(0))^/Gi=4™(t) 



Zs 



ET 



fc=-c 



ip{0)=0 



D[ip] cxp \ --S'setM (■ cos[ip{t) ~ ip{0)]. 



(112) 



This result may be used as a starting point for analyt- 
ical work and/or numerical calculations pl| ] . For further 
analysis, here we consider the semiclassical approxima- 
tion. 



B. Semiclassical limit 



^(-"+1)5^ ^,J^ 



h(3 



dada' a{a — a') 



(2m+l)!yo 
sinK(a-a')][C(a)-C(a')]'"+' (117) 



for odd orders and 

SET y (2to)! 7o ^ ' 



(118) 



For given winding number k, the path integral may be 
evaluated approximately by expanding around the clas- 
sical paths (^'^'^^(r) = if^ + Vj^T. An arbitrary path of 
winding number k may be written (/^(r) = (p^^'{t) + C,{t) 
with C(0) = C,(h(3) — 0. In terms of the Fourier coeffi- 
cients C,{vn) the action reads 

oo 

5set[^^'^ + C] = 2T,ikhn, + ^^e^tIC] + E '^'"^etK], 



(113) 



where the first term on the rhs is the topological contri- 
bution and 



"^set — "- 



7r2fc2 
pEc 



K 



the classical action of winding number fc, with the dimen- 
sionless parallel conductance g = G\\/Gk- The second 
order variational action 



'5^'S'sET = ^ XI ^SEtC*^") C(l^n) 



n=l 



(115) 



for even orders, with m 



1,2, 



Since A, 



setI^") is 



large for small (3Ec, the expansion (112) about the clas- 
sical path converges rapidly for high temperatures. At 
low temperatures Agg^p(^'„) vanishes for n < |fc| and the 
simple semiclassical approximation breaks down. The 
zero modes can be treated systematically for large g 
by considering quasi-classical trajectories with collective 
coordinates (sluggons) and fluctuations around them. 
This treatment, presented elsewhere [ p6| for the parti- 
tion function of the single electron box, lies outside the 
scope of the present work and we proceed with the high 
temperature expansion. Rewriting the cosine function in 
Eq. (112) as a sum of exponentials, we may perform the 
path integral and get for the correlator 



(114) (/i(t)/i(0))^/Gi 



fc— — OO 



27rikno 



exp 



-2E 



1 - COSJl^rnT) 
= 1 '^SET\'^in) 



i-^s: 



(k) 



(119) 



where the coefficients Ck read M 



is diagonal with the eigenvalues 
27r2n2 






pEc 



+ ge(n-\k\){n-\k\). 



The higher order terms in (|113|) read 



(116) 



Ck 



r(i 



fc+)r(i + fc_) 



r2(i-hA:)r(i 



(120) 



with fc± = fc + I ± i\/4Mfc + w2 and u = gf3Ec/2TT^. The 
contribution of the third order variational action cancels, 
thus the dominant correction to the semiclassical approx- 
imation stems from the fourth order term 



•^4 - T^gnP 2^ -(^j— (fe) ■ 



a.i^O 



(121) 



The corresponding expansion of the partition function 
^SET reads 



^SET — 



CSO 

\ '^ i^~f ,^2-Kiknq 



k— — oo 



\k\e 






(122) 
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with the same correction (121). The expansions ( |119| ) 
and ( |122| ) proceed in powers of (3Ec, however, terms in- 
volving u = gPEc/^Tr"^ are kept to all orders. This en- 
sures a meaningful result in the limit of moderately high 
temperatures also for large parallel conductance g. 



From Eq. (119) one obtains for the Fourier coefficients 



E{Vn 



hp 



dre"'"^(/i(T)/i(0)>^/Gi 



(123) 



the high temperature expansion 

47r " 



E{iy„ 



^SET 



Eq 



exp 



k— — oo 

oo 



k\e 



1 



-27Tiknn 



1=1 ^SEtI'^'7 



aivn+k) 



'n-\-k-\-m ) 






Ij: 



a{v, 



'n-t-fc+m+Z J 






1 



;.(fc) 



-5(z.„+fc)5r+0(/3^c)' 



(124) 



Since no 1/Aggrp(f„) term appears, the order of the ex- 
pression remains the same after analytical continuation. 
When E{vn) is analytically continued in the complex v 
plane, E{v) is analytic on each half plane Rei/>0 with 
a cut along the imaginary axis |4q ]. The representation 
of E{v) as a sum over winding numbers k shifts this cut 
to Re zy = fc for the fc'th term of the sum. Thus, in the 
phase representation, only the full sum shows the ana- 
lytic properties underlying the conductance formula ( 



The summands of the high temperature expansion (124) 
for winding number k are of the form g{k)f{\n ~\- k\) — 
e(-n - k)g{k)f{-n - k) + e{n + k)g{k)f{n + k), where 
/(n), g{n) are analytic functions and g{k) = g{—k). Now, 
the sum over winding numbers may be expressed as a sum 
over charges m 

J2 9{k)f{\n + k\) 

^ / dKe'^^"^^g{K)f{-n-K) 



fc — — OO 



+ dne''"'""'g{K)f{n + n). (125) 

Performing the limits ivn -^ lo -\- iS and w — > for the 
rhs of this equation and rewriting the result again as a 
sum over winding numbers, we obtain 



1 °° 

^™ni~I™ 1™ . E 9{k)f{\n + k\ 



fc— — c 



2^ 



E 9ik)j^fi\k + S\) 



dk 



'2^ 



.g(0)/'(0). (126) 



The sums in Eq. (124) may be performed exactly with 
the help of integral representations of the psi function 
and its derivatives fB2| 



%p{z) ^ I dt 



t 1 -e- 



and 



^(")(z) = (-1) 



71+1 



cxD j-n^ — zt 



dt- 



re- 



1-e- 



(127) 



(128) 



This way the high temperature expansion of the conduc- 
tance may be evaluated to read 

GsET^G,iZ-^^exp{-2[-f + i>{l + u)]/g} 

+ [ga{u) + t{u)] {pEcl2n^f+ 0{PEcf}. (129) 

The dependence on u = gf3Ec/2Tr^ is given in terms of 
two auxiliary functions 



a{u) 



7 -I- ■0(1 +-") -uil;'{l + u) 



{l-v)u 



and 



'(") = - 



377/.'(l + u)+ ^P{1 + u)[^ + 2z/''(l + u)] 



mi + u) + -ff 



+—iP{l + u)+ / dvE{u,v), 

with Lerch's transcendent (j){z,u,v) [Q and 
E{u,v) 



(131) 



1 



uln(w) (1 — v)u 



2v(j){v,2,l + u) 



1 — 2u" r 1 1 

ln(w)0(-,l,l + u)-f 0(-,2,l + u) 

V V V 



21n(u) ln(l -v) + - \n(v) + 3Li2(l - v) 



2(1 - w") 



[ln(l -v)+ v(t){v, 1, 1 + u)] }. (132) 



The high temperature expansion of Z is straightforward 
and reads 

ZsET = 1 + ga{u){l3Ec /27rY + 0{l3Ecf 

+2Ci cos(27rng) [l -f- 0{l3Ecf] , (133) 



which combines with Eq. (U29) to yield an analytical ex- 
pression for the high temperature conduction of a SET 
valid for arbitrary tunneling conductance. 
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C. Discussion of Results and Comparison with 
Experimental Data 

In Fig. |l^ the normalized conductance Gset/GcI is de- 
picted in dependence on the dimensionless gate voltage 
Ug for various temperatures f3Ec- The quantum correc- 
tions are more pronounced for lower temperatures where 
the gate voltage dependence becomes more significant. 
The oscillatory behavior of the conductance may be char- 
acterized in terms of a maximum Gmax = GsetU =1/2 
and minimum Gmin = GsetU =0 linear conductance. 



This gives for the conductance of a weakly conducting 
SET 



GSET 

Gel 



1 



l3Ec ^fl ^ 3C(3)\ 
x{PEcf+0{l3Ecf 



(135) 



in accordance with earlier work |55| - |57]| . In the region 
of strong tunneling, the quantity u is typically large even 
for the highest temperatures explored experimentally and 
the full expression (|129|), (|133|) must be used. 



g = 7.3 



0.64 



o 



O 0.6 



0.56 




/^c 




3 


^ 


3.5 


E 



4 





4.5 


E 



5 




5.5 





0.75 




FIG. 13. Conductance of a symmetrical SET in dependence 
on the dimensionless gate voltage rig for various temperatures 
and parallel conductance g = 7.3. 

We have compared our findings for the maximum and 
minimum as a function of temperature with recent ex- 
perimental data by Joyez et al. pq | for transistors with 
g — 0.6, 2.5 and 7.3. As seen from Fig. nj the theory 
describes the high temperature behavior of all junctions 
(results for g — 0.6 are not shown) down to tempera- 
tures where the current starts to modulate with the gate 
voltage. The parameters have not been adjusted to im- 
prove the fit but coincide with the values given in [M. 
The small deviations between theory and experiment for 
g = 7.3 near (3Ec = 1 may arise from experimental un- 
certainties in (3Ec [g^. We mention that the tempera- 
ture dependence of the conductance of the highly con- 
ducting SET {g = 7.3) is not within reach of previous 
theoretical work. The results obtained should be useful 
for experimental studies of even larger tunneling conduc- 
tances since the predictions remain valid for arbitrary 
values of g. 

In the region of weak tunneling, g < 1, the quantity u 
becomes small at high temperatures and we may replace 
a(u) and t{u) by 



FIG. 14. Maximum and minimum linear conductance in 
dependence on dimensionless temperature for two dimension- 
less parallel conductances g = 2.5 and 7.3 compared with 
experimental data (symbols) by Joyez et al. |l6|. 



VI. CONCLUSIONS 

In this article we have studied the conductance of 
nanofabricated metallic circuits showing Coulomb block- 
ade phenomena. We have treated electron tunneling non- 
perturbatively based on a path integral expression de- 
rived in Sec. II. Then, the frequency dependent linear 
conductance of a single tunnel junction embedded in an 
electromagnetic environment was calculated in the semi- 
classical approximation. We have shown that this ap- 
proximation is not only adequate for high temperatures 
but also in the limit of large conductance. As far as the 
leading quantum corrections are concerned, the tunnel 
junction was shown to be described as an effective linear 
element with an admittance that depends on the whole 
circuit. The predictions for the dc conductance were com- 
pared with recent experimental findings by two groups 
]l7| , p8[ and we found good agreement in the semiclassical 
regime of large conductance and/or high temperatures. 



cr(0) = 6C(3), t(0) = ttVIO. 



(134) 
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Further, we have shown that the low frequency behav- 
ior of the ac conductance can be calculated in terms of 
a renormalized capacitance which shows a linear depen- 
dence on the inverse temperature. 

In Sec. IV we applied the method to a linear array 
of N tunnel junctions and determined the effect of the 
environmental impedance on the conductance as well as 
the influence of the array length N. For large N the 
conductance dip becomes independent of the electromag- 
netic environment in accordance with previous calcula- 
tions p9| . For multi-junction circuits the configuration 
space of the phase variables was shown to be a torus, and 
contributions of nonvanishing winding numbers become 
relevant if one goes beyond the leading order quantum 
corrections. 

The conductance of the single electron transistor was 
determined in Sec. V by including nontrivial winding 
numbers leading to the gate voltage dependence of the 
conductance. The results were found to explain recent 
experimental data |16| for moderately low temperatures. 
For lower temperatures the contribution of sluggon tra- 
jectories [ p6[ has to be taken into account which was not 
elaborated here. 

The semiclassical theory presented has features in com- 
mon with the quasiclassical Langevin equation put for- 
ward in Ref. |3q]. The Gaussian approximation under- 
lying this approach is consistent with the semiclassical 
theory up to first order in /3Ec- Our results consistently 
include higher order terms in PEq for arbitrary tunnel- 
ing conductances. Non-Gaussian fluctuations are partic- 
ularly relevant in the moderately large tunneling regime. 
The analytical theory presented covers one edge in the 
temperature/conductance plane, arbitrary conductance 
and sufficiently high temperatures. Another edge is de- 
scribed by the perturbative approach, arbitrary temper- 
ature and sufficiently small tunneling conductance. Both 
theories arise naturally from the formally exact represen- 
tation of the current-current correlator which may also 
serve as a basis for Monte-Carlo simulations ||5^ that 
bridge between the semiclassical and perturbative results. 
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